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ABSTRACT 


On August 24 (UT) the Palomar Transient Factory (PTF) discovered PTFllkly 
(SN 201 lfe), the youngest and most nearby type la supernova (SN la) in decades. We 
followed this event up in the radio (centimeter and millimeter bands) and X-ray bands, 
starting about a day after the estimated explosion time. We present our analysis of 
the radio and X-ray observations, yielding the tightest constraints yet placed on the 
pre-explosion mass-loss rate from the progenitor system of this supernova. We find a 
robust limit of M ^ 10“ 8 (w/100kms~~ 1 ) M© yr _1 from sensitive X-ray non-detections, as 
well as a similar limit from radio data, which depends, however, on assumptions about 
microphysical parameters. We discuss our results in the context of single-degenerate 
models for SXe la and find that our observations modestly disfavor symbiotic progenitor 
models involving a red giant donor, but cannot constrain systems accreting from main- 
sequence or sub-giant stars, including the popular supersoft channel. In view of the 
proximity of PTFllkly and the sensitivity of our prompt observations we would have 
to wait for a long time (decade or longer) in order to more meaningfully probe the 
circumstellar matter of la supernovae. 


1. Introduction 

Type la supernovae (SNe) have served as an exquisite probe of cosmography (Riess et al. 
1998, Perlmutter et al. 1999) . As a result of this role, this class has been studied in unprecedented 
depth and breadth. Nonetheless the progenitors of la supernovae remain enigmatic. According 
to common wisdom a type la SN is due to the thermonuclear explosion of a white dwarf with 
a mass approaching the Chandrasekhar limit (Hoyle A Fowler 1960). The progenitor question 
is then one of understanding how a white dwarf can approach the Chandrasekhar mass (see 
review by Hillebrandt & Niemeyer 2000). In the single degenerate (SD) model (Whelan & Iben 
1973), a WD accretes mass from its hydrogen-rich star companion, reaches a mass close to the 
Chandrasekhar mass, which is sufficient to ignite carbon, and explodes. In the double degenerate 
(DD) model, a supernova results from the merger of two WDs (Iben A Tutukov 1984; Webbink 
1984; see also Yungelson & Livio 2000) . 
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It has long been suggested that radio and X-ray observations of type la SNe have the ability 
to provide diagnostics to distinguish between these two models (Boffi & Branch 1995, Eck et al. 
1995. Panagia et al. 2006). In most variations of the SD model, the winds from the donor star 
will enrich the circumstellar medium. The interaction of the blast wave from the supernova with 
the circumstellar medium can result in radio emission. In contrast, there is no expectation of 
circumstellar medium and hence of radio emission in the DD model. 

On UTC 2011 August 24.16 the Palomar Transient Factory (PTF; Law et al. 2009; Rau et al. 
2009) discovered PTFllklv, a rapidly rising transient, in the nearby (distance, d fa 6.4 Mpc; 
Shappee & Stanek 2011) galaxy Messier 101 (Nugent et al. 2011). Spectroscopy undertaken at the 
Liverpool Telescope led to a plausible la classification 1 and was soon confirmed by observations at 
the Lick 3-m telescope and the TNG. 

The apparent extra-ordinary youth and the proximity 2 of PTFllkly presents a unique 
opportunity to sensitivity probe the circumstellar medium of a type la supernova. Therefore, 
we immediately initiated (Gal-Yam et al. 2011) observations with the Swift Observatory, the 
Combined Array for Research in Millimeter-wave Astronomy (CARMA) 3 and the Expanded 
Very Large Array (EVLA) 4 . A few days later low frequency observations were undertaken at 
Westerbork Synthesis Radio Telescope (WSRT) 5 . X-ray observations were obtained with the Swift 

*at which point the event was rechristened to SN 201 lfe by the Central Bureau for Astronomical 
Telegrams 

2 The closest type la SN previous to PTF llkly is 1986G at a distance of 5.5 Mpc 

3 Support for CARMA construction was derived from the states of California, Illinois, and Mary- 
land. the James S. McDonnell Foundation, the Gordon and Betty Moore Foundation, the Kenneth 
T. and Eileen L. Norris Foundation, the University of Chicago, the Associates of the California 
Institute of Technology, and the National Science Foundation. Ongoing CARMA development and 
operations are supported by the National Science Foundation under a cooperative agreement, and 
by the CARMA partner universities. 

4 The EVLA is operated by the National Radio Astronomy Observatory, a facility of the National 
Science Foundation operated under cooperative agreement by Associated Universities. Inc. 

°The Westerbork Synthesis Radio Telescope is operated by ASTRON (Netherlands Foundation 
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and Chandra observatories. 


2. The Observations 

The early optical light curve of PTFllkly shows an extra-ordinary good fit to that expected 
from an exploding star (flux proportional to exponential of square of time). As a result the birth 
of the supernova can be accurately timed to a fraction an hour: UT 2011 August 23.69 (Nugent et 
al. 2011). Our first observations at both radio and X-ray bands were taken just over a day after 
the explosion. 

2.1. Radio 

The log of observations and the associated details can be found in Table 1. Our CARMA and 
EVLA observations include the earliest search for radio emission in cm-wave and mm-wave bands 
and subsequent observations include a very sensitive search in the 21-cm band obtained at the 
WSRT (see Figure 1). Following our first EVLA and CARMA observations, additional data was 
taken at a lower frequency (5 GHz) on UT 2011 August 25.8 by Chomiuk & Soderberg (2011). 
They reported a null detection with — 5 ± 6 /iJy. As can be gathered from Table 1 there are no 
detections at any epoch and in any band. In the next section we discuss the implications of these 
null detections. 


2.2. X-ray observation 

Following our classification of PTFllkly as a supernova we immediately triggered the Swift 
Observatory (see log of observations in Table 2) . The first Swift observations 6 began on UT 
2011 Aug 24.92. We used the Swift-X RT data products generator of the UK Swift Science Data 
Centre (SSDC: see Evans et al. 2009) to generate a single combined and astrometrically-corrected 

for Radio Astronomy) with support from the Netherlands Organization for Scientific Research 
N WO : . 

6 Target ID 32081. with initial ObsID 32081001 lasting 4.5 ksec over three consecutive orbits, 
followed by 53.7 ksec of observations through Sep 8.58 (ObsIDs 32081002 32081029). 



Table 1: Log of radio observations 


Start 

UT 

AT 

day 

r 

min 

Facility 

V 

GHz 

s„ 

+Jy 

Luminosity 
^ 10 24 ergs _1 Hz -1 

M 

10 V/-7 < .1 M q yr“ 

Note 

i 

Aug 24.98 

1.4 

178 

CARMA 

93 

-16 + 510 

75 

23 

(i) 

Aug 25.02 

1.3 

37 

EVLA 

8.5 

-4.4 + 25 

3.7 

1.7 

(2) 

Aug 27.71 

4.0 

35 

EVLA 

5.9 

1.3 + 7 

1.0 

1.1 

(3) 

Aug 28.25 

4.8 

800 

WSRT 

4.9 

-71 ± 34 

5.0 

2.6 

(4) 

Aug 29.97 

6.3 

29 

EVLA 

5.9 

- 0.9 + 9 

1.3 

3.1 

(3) 

Aug 31.38 

8 

630 

WSRT 

1.4 

2 + 25 

3.7 

1.7 

(4) 

Notes: The 

columns 

starting from left 

to right are as follows: start 

of integration in 

UT ; mean 


epoch of observation (in days since explosion); integration time in minutes; facility; central frequency 
in GHz: nominal flux and associated rms in the vicinity of PTFllkly in /rJy; the corresponding 3-<r 
spectral luminosity assuming a distance of 6.4 Mpc to M101; inferred upper limit to the mass loss rate 
(see §3.1 for explanation of parameters); specific notes (see below). Following Nugent et al. (2011) we 
assume that the explosion time of PTFllkly is UT 2011 August 23.69. The following packages were 
employed to reduce the data AIPS (EVLA), Miriad (CARMA) and NEWSTAR (WSRT). (1) The CARMA 
observations were obtained in the E configuration using only nine 6-m antennas. The larger antennas were 
not available owing to reconfiguration of the array. Bandwidth of 8000 MHz. Calibrators: J1 153+495 & 
J1642+689 (phase) and MWC349 (flux). (2) Bandwidth of 256 MHz. Calibrators: J1419+5423 (phase) 
and 3C286 (flux). (3) Band width of 2000 MHz. Data obtained under Director’s discretionary time (PI: 
A. Soderberg). (4) Band width of 8 x 20 MHz. For each of the eight IF channels, only 3/4 of the channel 
bandwidth was used in making the map. The flux density calibrators used were 3C147 and 3C286 on the 
Baars et al. (1977) scale. 
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event file from the first 58.2 ksec of exposure. Six events are found within an aperture of radius 9 
XRT pixels (21.2") centered on the position of PTFllkly. This can be compared to independent 
background expectations of 7.0 ± 0.4 and 11.7 ± 0.7 counts, respectively. The corresponding 
90%-confidence limits on total source counts and aperture-corrected average count rate are n 7 < 4.3 
and r x <0.1 x 10" 3 s" 1 , respectively. 

However, we note that two of the counts within the source region arrive within the first 2 ksec 
- an a ■priori unlikely occurrence (at ~90% level of confidence) - and could be taken as evidence of 
early, bright X-ray emission from PTF llkly, at the L x ~ 10' 18 ergs” 1 level (see also Fox 2011a, b). 
We realize that this statistical evidence does not warrant a claim of detection. Nonetheless, given 
that the observations were done at an extraordinarily early epoch of 1.21 day post-explosion we 
think it worth highlighting this issue for the benefit of future observers. 

Our upper limit of n 7 < 4.8 counts, from the first Swift observation sequence (4.5 ksec 
exposure) alone, provides the following upper limits on X-ray flux (over the energy range 
0.3-10 keV) 7 : F x < 6.2 x 10“ 14 erg cm -2 s _1 (thermal bremsstrahlung model with kT = lOkeV) 
and F x < 5.0 x 10~ 14 ergcm~ 2 s _1 (power law model with photon index, F = 2). Corresponding 
upper limits to the X-ray luminosity are L x < 3.0 x 10 38 erg s~~ l and L x < 2.5 x 10 38 ergs _1 , 
respectively. 

Observations 8 with the Chandra X-ray Observatory began on UT 2011 Aug 27.44 and 
lasted for 49.7ksec; the mean epoch is UT 2011 Aug 27.74 (corresponding to 4days post- 
explosion) (Hughes et al. 2011). No soft proton flaring was evident during the observation. Only 
one event was found in the source region, a 2.5-arcsec aperture centered on the SN. whereas 
2.23 ±0.1 background counts were expected (in an equivalent aperture). The 90%-confidence upper 
limit on the expectancy value of a Poisson process that generates one count is 3.9 counts. Ignoring 
the background contribution, after applying the aperture correction (as prescribed by Feigelson et 
al. 2002). we find an upper limit to the count rate, r x ^0.11 x 10~ 3 s“ 1 (90%-confidence). 

‘ We assume a Galactic ISM contribution of iV# = 1.8 x 10 2O cm” 2 ; Kalberla et al. 2005. There 
may be an additional similar contribution from M101, depending on the depth of the line of sight 
to the SN; Kamphuis et al. (1991). 

8 ObsID 14341; PI J. Hughes 
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For the thermal bremsstrahlung model this upper limit translates to F x < 9.0 x 
10 -16 erg cm -2 s -1 (0.3 -8.0 keV), corresponding to L x < 4.4 x 10 36 ergs -1 . For the power law 
model we find F x < 8.2 x 10 -16 erg cm -2 s -1 (0. 3-8.0 keV). corresponding to L x < 4.0 x 10 36 ergs -1 . 

3. Mass-loss rate constraints 

The theory of radio emission from SNe relevant to Type I events is discussed in Chevalier 
(1982, 1998) and Chevalier & Fransson (2006) and a summary of radio observations of supernovae 
can be found in Weiler et al. (2002). Panagia et al. (2006) report a comprehensive summary of 
searches for radio emission from la supernovae. 

The basic physics is as follows. The supernova shock-wave ploughs through the circumstellar 
medium (CSM). In the post-shock layer, electrons are accelerated to relativistic speeds and strong 
magnetic fields also appear to be generated. The relativistic electrons then radiate in the radio 
via synchrotron emission. X-rays are emitted via two energy channels: thermal bremsstrahlung 
emission from hot post-shocked gas and inverse Compton scattering of the optical photons from 
the supernova by the relativistic electrons. 

The density of the circumstellar medium is a key physical parameter. After all a strong shock 
can only be generated if the SN blast wave can run into CSM. Thus, the strength of the radio and 
X-ray emission allows us to diagnose the CSM density. 

For a star which has been losing matter at a constant rate, M, the circumstellar density has 
the following radial profile: p(r) = n(r)p = M / (4irr 2 w) where w is the wind velocity, r is the 
radial distance from the star, n is the particle density and p is the mean atomic weight of the 
circumstellar matter. We assume that the circumstellar matter is ionized (by the shock breakout 
and by radiation from the young supernova). 

3.1. Constraints from Radio Measurements 

Massive stars exploding as type Ib/Ic SNe with their fast moving ejecta and low density CSM 
provide a good starting point to discuss radio emission from la supernovae. The spectrum of the 
radio emission from type Ib/Ic events follows the “Synchrotron Self-Absorption” (SSA) form: a 
v 5 ' 2 low-frequency tail (optically thick regime) and a declining power law, v a at high frequency 
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(optically thin regime). For most well studied SNe a ~ — 1. Diagnosis of the CSM is based on the 
peak frequency (u m : the synchrotron optical depth is unity at this frequency) and the peak flux 
(S m ) and the evolution thereof. 

We use the basic SSA formulation as in Chevalier (1998) with some modifications (noted 
below). We assume a fraction of electrons are accelerated to relativistic energies and with a power 
law distribution, dN/dE = NqE~ p ; here E = ' ym e c 2 is the energy of electrons and 7 is the Lorentz 
factor. Electrons which are not relativistic will not radiate strongly and thus we introduce a 
minimum Lorentz factor, 7 min . Only electrons with 7 > 7 m j n are assumed to contribute to the 
radio emission. 

As in all strong shocks, in addition to acceleration of electrons, strong magnetic fields 
appear to be generated in the post shock gas. Our current understanding is such that we simply 
parametrize the relative energy fractions of each of these two components. 

The thermal energy density of the post-shocked gas is | pv 2 (where v s is the shock speed). 
The magnetic energy density is B 2 /(8tt). We denote the ratio of this energy density to that of the 
thermal energy density of post-shocked gas by e B - We find 


u B 


B 2 

8n 


6.8 x 10 9 m 2 e s 


( MtMsyi- 1 ] 
\ to [km s" 1 ] 


t d 2 erg cm 3 ; 


( 1 ) 


here, m is the power law exponent in the equation relating the radius of the shock to time, R s oc t m 
and 0 is the time post-explosion in days. In a similar manner, we let e e denote the ratio of the 
energy density of relativistic electrons (7 > 7 min ) to the thermal energy density. 

The use of radio diagnostics (unfortunately) involves the values of several additional 
parameters. The first parameter is p, the power law index of the relativistic electrons. Theory 
and observations of this value agree that this value should be p « 3 (see Weiler et al. 2002; 
Chevalier & Fransson 2006). Next, we make the simplifying assumption that the blast wave moves 
at a constant velocity, v s ; this is equivalent to setting m — 1. This is admittedly a simplification 9 . 


9 We are interested in the very outer layers of the exploded star. Note that with (say) a mass 
loss rate of lO -8 M 0 yr -1 and a wind velocity of 10' ems” 1 the leading edge of the shock (v s = 
4 x 10 9 cm/s. say) would have swept up a mere 10~ 8 M© by day 1. It is the velocity-density structure 
of such thin outermost layer that determines v s . 
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However, given our non-detections we felt that an analysis more sophisticated than this was not 
warranted. 

Finally, we come to the microphysics parameters, e e and e B . It appears from well studied 
SN shocks that e e is usually about 0.1 with modest dispersion (Chevalier & Fransson 2006). 

In contrast, e B appears to be a highly variable parameter. This is quite understandable from 
simple considerations. The magnetic field is generated both in the post-shock gas as well as by 
compression of the field already present in the pre-shock gas. For instance a radio and X-ray 
modeling of SN 2002ap (a type Ic supernova with a presumed Wolf-Rayet progenitor) by Bjornsson 
& Fransson (2004) find that the assumption of equipartition (e B — e e ) lead to inferring M much 
lower than the expected for a Wolf-Rayet progenitor. The inferred mass loss rate can be increased if 
6 b = 2 x 10“ 3 . We find a similar imbalance between e B and e e for SN 201 ldh (a type lib supernova; 
Horesh et al. 2011; see also Soderberg et al. 2011). For this reason, unless stated otherwise, we 
will fix e e = 0.1 and assume that e B is a free parameter. Finally, since only relativistic electrons 
contribute to radio and X-ray emission we adopt 7 min = 2. 

Following Chevalier (1998), the radio spectral luminosity in the optically-thin regime is given 


by 


L„ = 1 " 


2c, 


-(p- 1)/2 


ergs ] Hz x , 


and the synchrotron-self-absorbed luminosity is 


L u = 4 tt 2 fR^B^l-f 

Cg \ — C] 


5/2 


ergs Hz 


( 2 ) 


( 3 ) 


where R s = v s t is the shock-wave radius, B is the strength of the magnetic field; and the constants 
Ci, c 5 , and eg can be found in Pacholczvk (1970). N 0 can be straightforwardly shown to be 
e e /(87r e B )B 2 (p — 2)(7 mm m e c 2 ) p ~ 2 . Equation 2 can be simplified to yield 


L y = 3 x 10' 


34 1287r 3 /u 3 




syr 


u; [km s" 1 ] 


v 

2ci 


fjeres J Hz . 


( 4 ) 


Given that the earliest observation(s) were undertaken only a day after the explosion it is 
prudent to check if there is significant free-free absorption. The free-free optical depth is 


r ff = 3.3 X 10~~ 7 rr L35 r / GH2 EM 


( 5 ) 
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where T = 10 4 T 4 is the electron temperature (in degrees Kelvin) 10 , z'ghz is the frequency in GHz 
and EM is the emission measure, the integral of n‘f. along the line of sight, and in units of cm -6 pc. 
The emission measure from a radius, say, r* to infinity is 



where n* is the density of electrons at radius r*. Putting these equations together the free- free 
optical depth is 

r ff » 0.5 (7) 

For PTFllkly, the photospheric velocity is v = 2 x 10 9 cms -1 . The blast wave will at 
least have this velocity and likely twice this value (Chevalier & Fransson 2006; Fryer et al. 2007; 
Soderberg et al. 2010). For representative values (M = 3 x 10~ 7 M e yr -1 , w 7 = 1, Vg = 4) we 
find the following optical depth at the first epoch of our observations: T ff [z/GHz = 1.4,5,8,95] = 
[3.5,0.24,0.09,0]. For interesting values of M (^ 10~ 8 M 0 yr _1 ) free-free optical depth is not 
important for the observations reported here. 

Equation 4 provides the starting point for the discussion. This equation shows that the 
spectral luminosity can constrain M/w provided that we have a good grasp of the blast wave 
dynamics and microphysics of particle acceleration and magnetic field generation. We have argued 
above that p ~ 3. 6 e ~ 0.1 and v s «4x 10 9 crns F We adopt these values and proceed with the 
discussion. 

To start with we can see that the observations reported here (Table 1) yield the lowest limits 
on the radio luminosity of very young la supernovae, L u £ 10 24 erg s _1 Hz _ F This limit then 
constrains the following parameter, B = eM/iu, where e = ^/e B e e . We have deliberately not quoted 
the limits on M from previous literature since radio measurements yield not M but B and this 
quantity depends on the unknown parameter, e# which (as summarized above) can vary by orders 
of magnitude. Despite this it is clear that the observations reported here present the most sensitive 
limits to M to date (Figure 2). 

From Figure 2 we note that B can be constrained as follo-ws: the optically thin regime offers a 
10 We use the convention of X n = X/I0 n where it is assumed, unless explicitly specified, that the 


units are CGS. 
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lower bound whereas the optically thick regime (SSA and free-free) an upper bound. The upper 
bound is not interesting since the simplest explanation for the absence of radio emission is that the 
explosion takes place in a vacuum (or very low density circumstellar medium). Besides, from past 
studies, there is no indication of M in the range of 10" 6 M® yr" 1 that is indicated from the upper 
bounds. 

The constraint deduced for each observation are summarized in Table 1. Combining all the 
constraints we find M ^lx 10~ 8 u>7(0.1/e)M©yr^ 1 . 

3.2. Constraints from X-ray data 

Past X-ray observations, typically undertaken no earlier than a week past the explosion, have 
resulted in upper limits at the level of Lx 10 39 ergs" 1 . A claim of detection of emission from 
SN2005ke (Lx ~ 2 x 10 38 ergs _1 ; Immler et al. 2006) has been disputed by Hughes et al. (2007). 
Immler & Russell (2011) reported a detection in a 2-month stack of Swift XRT data of SN 2011by. 
However, subsequent high angular resolution Chandra observations by Pooley (2011) strongly 
suggest that the emission arose from a steady source 2.3-arcsecond away from the supernova. The 
most sensitive and constraining observation (prior to this event) was from a 20-ksec Chandra 
observation of SN 2002bo, obtained 9 days past the explosion (Hughes et al. 2007). The 3-cr upper 
limit was L x < 1-3 x 10 38 ergs _1 (2-10 keV). 

The Swift and Chandra X-ray observations of PTF llkly provide the earliest (1 day) as well 
as the most sensitive limits («30 times deeper relative to previous upper limits) on X-ray emission 
from SNe la. As noted earlier (§3) there are two durable sources of X-ray emission: thermal 
bremsstrahlung and inverse Compton (IC) scattering of the SN photons by relativistic electrons. 
For such strong shocks, kT=10keV is a characteristic temperature (bremsstrahlung). IC scattering 
from 1 eV to the 0.1-10 keV band requires electrons with Lorentz factor 7 in the range of 10 to 
100, a reasonable range for strong (non-relativistic) shocks. 

We first consider a model in which we ascribe all X-ray emission (if any) to ther- 
mal bremsstrahlung. Following the model discussed in Immler et al. (2006), Lx = 

(Trm 2 )" 1 A(T) (M/w 2 (vt)~ l , and with the assumptions of an H+He plasma as in §2. 

^reverse = 30 Lforward; and A (T) = 3x 10" 23 erg cm 3 s _1 , the upper limit of L© < 4.4 x lO'^ergs" 1 
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Fig. 2. — The model light curves for four different values of B = eM/w where e 2 — eBe e , w is 

the wind velocity and M is the mass loss rate (see also §3.1). We assume the following normal- 
ization. e — 0.1. w — 10'cms -1 and v s — 4 x 10 9 cms“ 1 . Model light curves for four values of 
B = [10~ 6 , 10" ' . 10~ 8 . 10 _9 ](0.1/e)tO7 M^yr -1 are shown. Top panel is for X band (8 GHz) while bot- 
tom panel is for C (5 GHz) band. The open triangles in both panels shov r upper limits from Panagia et al. 
(2006). The gray triangle in the top panel is the best limit from Panagia et al. (2006) while that in the 
bottom panel is the upper limit from Hancock et al, (2011). Red triangles are measurements of PTFllkly 
and presented in Table 1 . 
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translates into a progenitor mass-loss constraint of M ^5 x 10~ 6 W7 yr -1 . Scaling from 

the upper limit for SN 2002bo, given the distinct and more conservative model assumptions of 
Hughes et al. (2007), suggests a limit for PTF llkly of M ^4x 10“ 5 w 7 Af 0 yr _1 . 

Next, we consider the IC scattering model. From Chevalier & Fransson (2006) we find that 
the IC luminosity is 

FlC « 10 36 7min(^j)v9^ 1 M-8^f 1 ( 10 42^g S -l ) ^ S ~ 1 ’ ( 8 ) 

where Lsn is the bolometric luminosity from the supernova. Clearly, the IC luminosity is 
independent of the highly variable e B and, unlike the radio luminosity (Equation 4) has a gentle 
dependence on v s . Assuming e e — 0.1 and noting that L s n ~ 2 x 10 42 ergs“ 1 on day 4 we find 
M < 0 10 _8 tn7 Mq yr _1 (90% confidence limit). Unlike the radio case this is a robust limit since it 
does not depend on eg. 


4. Discussion 

As noted in §1 the commonly accepted SN la model is a thermonuclear explosion of a 
Carbon-h Oxygen white dwarf close to the Chandrasekhar mass. However, such white dwarfs have 
masses smaller than the Chandrasekhar mass and so the la models necessarily involve the growth 
of mass of the white dwarf. One such model is the coalescence of tw r o white dwarfs whose total 
mass is in the vicinity of the Chandrasekhar mass (see van Kerkwijk et al. 2010 and references 
therein). There is no expectation of an enriched circumstellar medium for (long lived) double 
degenerate (DD) systems. Thus the blast wave from the supernova is not expected to produce any 
strong emission either in radio or X-ray bands. 

On the other hand, in the “Single Degenerate” (SD) model the white dwarf can grow by a 
steady transfer of matter from a companion. This can occur from a close binary in Roche Lobe 
Overflow (RLOF), or by accretion directly from the wind of a companion (the symbiotic channel). 
In this class of models one may expect a circumstellar medium enriched by matter from the donor 
star either directly (by a wind from the donor star) or indirectly (matter that the white dwarf is 
unable to accrete is expelled from the binary system). 

It has long been the hope that early and sensitive radio and X-ray observations of nearby la 
supernovae would diagnose the close in circumstellar matter and thereby discriminate between the 
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SD and DD models. 

Here, we report the earliest radio and X-ray observations of PTFllkly and sensitive radio and 
X-ray observations undertaken by us and others at premier radio and X-ray facilities. Despite the 
sensitivity and rapidity no signal is seen at either radio or X-ray bands. 

The absence of the signal is consistent with the expectation of the DD model but would hardly 
constitute proof of this channel. Below, we confront popular SD models with our null detection. 

The radio observations find that M <; 10~ 8 u;7(0.1/e) M©yr -1 , where w = 10' uycms -1 is 
the velocity with which the matter is ejected from the binary system and e = yTieJ. The X-ray 
observations yield an upper limit of 10~ 8 wf 1 (0.1/e e ) M© yr _1 . To our knowledge, these are the 
most sensitive limits on M reported for any la supernova, to date. 

Returning to the SD model, the growth in the mass of the white dwarf is not assured. To start 
with, the mass transfer rate has to be high enough to prevent a nova explosion carrying off all 
the accreted matter. At high enough rates, steady burning can occur on the white dwarf and the 
white dwarf can grow in mass. Accretion at a rate of rs 10“' M©yr -1 is thought to be necessary 
for stable accretion and nuclear burning on the surface of a white dwarf (Nomoto 1982). 

Systems in RLOF may achieve the required accretion rate - the class of binaries known as 
’supersoft x-ray sources’ (van de Heuvel et al. 1992 ) show such behavior. But in addition to steady 
accreation, mass loss from the system is also thought to occur. Figure 1 from Han & Podsiadlowski 
(2004) gives theoretically expected mass-loss rates for the classical supersoft channel using the WD 
accretion efficiencies from Hachisu et al. (1999). In this phase, the mass loss from the system can 
reach 10“ 6 M©yr“ 1 and, in extreme (but rare) cases, more than 10 _5 M©yr _1 . The velocity of this 
material depends on the details of the ejection process; if it comes from the neighborhood of the 
WD, it is expected to be several 10 8 cms _1 ; if it comes from a circum-binary envelope, its velocity 
should be at least of order the typical orbital velocity, i.e. a few 10 7 cms _1 . Note, however, that in 
the vast majority of cases, the binary will not be in this wind phase at the time of explosion, and 
only a fraction of the mass transferred may be lost from the system at that time; unfortunately, 
the exact amount is not well constrained by the theoretical models (perhaps a few 10 ~ 8 M© with a 
velocity of at least a few 10 ‘ cms -1 ). 

The radio limits for PTF llkly do not rule out mass loss at that level. If the donor star is a 
main-sequence star or a slightly evolved sub-giant, its stellar wind is unlikely to be an important 
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contribution to the systemic mass loss. 

In the ‘symbiotic’ SN la channel, the stellar wind from the donor star may be more important 
than mass loss associated with the mass-transfer process. If the companion is a Mira variable (in a 
D-type symbiotic), the expected mass loss is 10" 6 to 10“ 4 M©yr _1 (e.g., Zijlstra 2006) with very 
low terminal velocity (w ~ 10 km s'" 1 or even less). Such high mass loss can clearly be ruled out in 
the case of PTF llkly. 

In an S-type symbiotic, which is a more likely SN la progenitor (Hachisu et al. 1999), the 
donor star is a much less evolved red giant with much lower mass loss. Scaquist & Taylor 
(1990) estimate typical rates of 10~ 8 -10~ 7 M 0 yr -1 with large modeling uncertainties. A case of 
particular interest is the symbiotic binary EG And which contains a M2. 4 red giant that is very 
similar to the red giant in the SN la candidate system RS Oph. The wind from EG And has a 
measured high terminal velocity of 75 km s 1 (Espey & Crowley 2008) with an uncertain estimate 
of the mass loss rate of about 10 _8 M o yr _1 (Crowley 2006). For comparison, the observed wind 
features in RS Oph have velocities in the range of 30-60 km s _1 (Patat et al. 2011). In addition, 
observations of RS Oph show that the mass loss is very asymmetric and is strongly confined to 
the orbital plane (O’Brien et al. 2006), consistent with hvdrodynamical modeling of the mass 
loss from such systems (Walder et al. 2008; Mohamed, Booth & Podsiadlowski [in preparation]). 
This introduces a possibly important viewing dependence of the radio signal. In this context, it 
is worth noting that these systems (both in the supersoft and the symbiotic channel) experience 
recurrent novae, typically every 10-20 yr when the white dwarf is close to the Chandrasekhar mass 
(Schaefer 2010), which may produce low-density cavities in the immediate circumstellar medium 
(Wood-Vasey & Sokoloski 2006). The upper limit for the mass- loss rate deduced from our radio 
limits for the progenitor of PTF llkly is comparable to the expected mass loss from a symbiotic 
like RS Oph. However, considering the uncertainties noted above, it may be premature to rule out 
such a system on the basis of these observations alone. 

The limits we have derived from radio and X-ray observations on the progenitor system 
of PTF llkly can be placed into context with complementary constraints from a variety of 
independent techniques. Li et al. (2011) examine pre-explosion HST images of M101 and derive 
limits on the brightness of any mass-donating companion. The lack of any optical emission at the 
location of PTF l lkly directly rules out luminous red giants and the vast majority of He stars. 
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consistent with the limits we have derived here. 

In summary, we present the most sensitive early radio and X-ray observations to date. Despite 
the rapid response and excellent sensitivity our observations can constrain only the symbiotic SD 
channel for PTFllkly, but cannot rule out a main-sequence or sub-giant donor channel. 

Prior observations which were both less sensitive and at later times nominally came to the 
same conclusion. However, the inference of M from radio data depends on two microphysics 
parameters, e e , (b and two velocities, v s (shock speed) and w, the velocity with which matter is 
ejected from the binary system. In principle, v s can be estimated from photospheric velocities and 
a model for the exploding white dwarf. The value of e e empirically shows modest variation and one 
can assume e e « 0.1. However, e B , on empirical and theoretical grounds can vary tremendously 
(with €b ~ 0.1 being a maximum plausible value). 

In the past, particularly in the radio literature, the assumed values were (in our opinion) 
rather optimistic: w ~ 10 6 cms -1 and e B = 0.1. Using the maximum possible value for e B ~ 0.1 we 
find from our radio observations that M <; 10 Ur- M ( . yr _l . Smaller values of e B will only make 
this limit worse (proportionally larger). Given that PTFllkly was one of the closest la supernovae 
it is not likely that the limits presented here would be bested in the near term. 

The X-ray observations provide a less model dependent (and hence more robust) estimate 
of M. The X-ray observations (especially if undertaken at peak) tightly constrain Me e /w. The 
microphysical parameter, e e has far less dispersion as compared to e B and as such the X-ray 
observations yield a robust estimate of M (relative to that obtained from radio observations). We 
find M 10 -8 u;7 M 0 yr -1 . The X-ray observations are quite sensitive and it is not likely that 
these limits will be easily surpassed in this decade. 
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Table 2: Log of X-ray observations 


Start 

AT 

T 

Facility 

Band 

F x 

Luminosity 

M 

UT 

day 

ksec 


keV 

g 10" lc ergc??r -2 s" 1 

10 36 ergs _1 

10~ H wr/e e -i Mq yr“ J 

Aug 24.92 

1.21 

4.5 

Swift 

[0.3-10] 

500 

250 

20 

Aug 27.44 

4 

49.7 

Chandra 

[0.3-8] 

8.2 

4 

1.1 


Notes: The columns starting from left to right are as follows: start of integration in UT; mean 

epoch of observation (in days since explosion); integration time in minutes; facility; Energy band in keV; 
flux in 10 _16 erg cm~ 2 s -1 ; the corresponding luminosity limit assuming a power law model with photon 
index T = 2 and a distance of 6.4 Mpc to M101; inferred upper limit to the mass loss rate (see §3.2 for 
explanation of parameters). 
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